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ABSTRACT 



Context. According to the current cosmological paradigm, large scale structures form hierarchically in the Universe. Clusters of 
galaxies grow through a continuous accretion of mass, and the presence of cluster substructures can be interpreted as signature of this 
process. Nevertheless, the rate and manner of mass accretion events are still matters of debate. 

Aims. We have analysed the presence of substructures in one of the largest sample of nearby cluster galaxies available in the literature. 
We have determined the fraction of clusters with substructure and the properties of the galaxies located in such substructures. 
Methods. Substructure in the galaxy clusters was studied using the Dressler-Shectman test, which was calibrated through extensive 
Monte Carlo simulations of galaxy clusters similar to real ones. In order to avoid possible biases in the results due to differing 
incompleteness among clusters, we selected two galaxy populations: a) galaxies brighter than M^ =-20 located in clusters at z < 0. 1 
(ECl); and b) galaxies of brightness M, < -19 located at z < 0.07 (EC2). 

Results. In the inner cluster regions (r < r2oo) 11% and 33% of the clusters of ECl and EC2 respectively show substructure. This 
fraction is larger in the outer cluster regions (a; 55%) for ECl and EC2 samples. Cluster global properties, such as cTc, ft or Ami2 
do not depend on the amount of cluster substructure. We have studied the properties of individual galaxies located in substructures 
in the ECl and EC2 galaxy populations. The fraction of galaxies within substructures is larger in the outer cluster regions when 
fainter galaxies are included. The distribution of relative velocities of galaxies within substructures suggest that they consist of an 
infalling population mixed with backsplash galaxies. We can not rule out that the infalling galaxy population located in substructures 
are genuine field ones. 
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1. Introduction 

Pioneering photographic surveys (e.g. lAbell 1958t 
IZwickv 19681 : IShectman 1985h showed that galaxies were 



distributed not uniformly but in clumps, forming aggregations 
called galaxy clusters. These gravitationally bound systems 
are (together with galaxy superclusters) the largest and most 
massive structures known in the Universe. These large galaxy 
aggr egations have been fo rmed by merging smaller galaxy units 
(e.g. lWhite & Rees 1978h . 

The increasing availability of galaxy surveys during 
the last decade has shown that a considerable fraction of 
nearby galaxy clust ers are not yet virialized systems (see e.g. 
iRamella et al. 20071 and references therein). The most obvious 
observational signature is the presence of multi-modality in 
the spatial and velocity distributions of galaxies and/or gas. 
Substructure can be defined as the presence of two or more 
galaxy aggregations within a cluster of galaxies. These galaxy 
aggregations could be associated with late cluster formation 
due to a recent merger of small groups forming a cluster, 
secondary infalls of galaxy groups in a virialized cluster, or 
gravitationally bound subclusters still in an infalling process. 
Studying cluster substructure therefore allows us to investigate 
the process by which clusters form. Indeed, substructure in 
galaxy clusters depends on the expansion rate of the Universe, 
and its study can also constrain the cosmological model for 
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Structure for mation an d our understanding of the nature of 
dark matter dRichstone et al. 1992; Kauffmann & White 1993; 
'Lacev & Col e 1 9931; iMohr et al . 1995; Thomas et al. 1998 ; 
Markevitch et al. 2004 IClowe et al. 200^. Moreover, the 
identification of cluster substructure is an important is- 
sue, since they can strongly affect estimates of global 
cluster parameters such as mass and velocity dispe rsion 
dPereaetal. 1990HEscaleraetal. 1994l;lGirardietal. 19961) . 

The presence of substructure in clusters can also influence 
the evolution of the constituent galaxies. Galaxy clusters 
are popular sites for studying environmental influences on 
galaxy evolution. Cluster galaxies are substantially different 
from galaxies in the field. This idea is supported by many 
observational results, such as the morphology-density relation 
(Dressier 198 3), the low fraction of star- forming galaxies 
(e.g. Popess o et al. 2007h . or the different structural param- 
eters shown by gal axies in hig h density en vironments with 
respect to the field ( Aguerri et a l. 2004; Gut ierrez et al. 2004 ; 
Barazza et al. 2009'; 'Marinova et al. 2009; Aguerri et al. 20091: 
Aguerri & Gonzalez-Garcia 2009). It was previously thought 
that cluster and field galaxies had different origins. Nowadays, 
these differences are explained in terms of the different evo- 
lution of the high density and the field galaxy population. 
Hence, all previous observational results could be explained 
in terms of several large- and small-scale mechanisms present 
in high density environments and not i n the field, amon g 
which we can mention galaxy harassment dMoore et al. 19961) . 
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ram-pressure st ripping ( Gunn & Gott 1972t lOuilis et al. 2000h . 
tidal stripping (Merrittl98J; 'Merritt 1984"), major and mino r 
galaxy mergers ( Aguerri et al. 2001 ; Eliche-Moral et al. 2006h . 
and starvation ( Larson etal. 19801) . Within this framework 
of galaxy evolution, where the environemnt plays an im- 
portant role, substructure should be taken into account in 
seeking a global picture of the evolution of galaxies. There 
is thus observational evidence that a higher fraction of clus- 
ter galaxies with ongoing starburst episodes are located 
in substructures or in the region s of cluster-subcluster 



interactions dC aldwel let al. 19931: iBiviano etal. 19971: 



iMoss & Whittle 20bot IPoggianti et al. 20 04^ Moreover, it 
has been observed that the presence of substructures is also 
directly related to the growth of the brightest cluster galaxies 
(BCGs). Thus, brighte r BCGs ar e located in clusters with less 
substructure (Ramella et al. 2007h . 

The analysis of substructure can be performed using 
the projected phase-s pace distribution of cluster galaxies 
(e.g. iDressler & Shectman 1988) . the surface-brightness dis- 
tribut ion and tempera ture of the X-ray emitting intra-cluster 
gas dBriel etal. 19 92*). or the shear pattern in the back- 
ground galaxy dis tribution induced by gravitational lensing 
dAbdelsalam et al. 1998 ). In this paper we focus on the study 
of the phase-space distribution of galaxies in clusters, tradi- 
tionally studied using different statistical tests. Depending on 
the information used, those statistical tests can be classified as 
ID, 2D or 3D. The ID tests use only the velocity information. 
The null hypothesis is that the velocity distribution in relaxed 
galaxy clusters should be Gaussian. Thus, these tests study the 
degree of non-gaussianity of the velocity distribution of galaxies 
analysing (among other things) its skewness, kurtosis, asymme- 
try index or tidal index (Beers et al. 1990). The 2D tests assume 
that the galaxy distribution in virialized clusters should follow a 
Poissonian distribution. The most popular 2D tests are, among 
others: the an gular separation test dWest et al. 1995), the ;6-test 
dWestetal. 19 88). the Fourier elongation test (Mohr et al. 1993") 
and the Lee statistic (Fitchett & Webster 1987; Fitchett 1988). 
The 3D tests use the spatial plus the radial velocity infor- 
mation. The null hypothesis of these tests is that no corre- 
lation between position and velocity exists. The most popu- 
lar of these tests are the Lee 3D test, the A-test and the a- 
test (IWest & Bothun 19901) . Numerical simulations have shown 
that 3D tests are the most efficient for det ecting substructure 
dKnebe & Miiller 20001: IPinknev et al. 19961) . 

The previously mentioned tests have been applied to 
different cluster samples and have revealed that about 1/3 
of the nearby g alaxy clust ers are not yet v iri ahzed sys- 
tems (iWest & Bothun 1 990t iRhee et al. 1 99 ll: iBird 1 994 : 



Escalera etal. 1994 IWestetal. 19951: iGirardi etal. 1997 
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Burget t et al. 2004 

Most 



Solanes et al. 1999'; 'Biviano et al. 2002"; 
Flin & Krywult 2006; Ramella et al. 2007). Most of these 
studies were done on small galaxy cluster samples and/or they 
contain no velocity information on the galaxies. Few nearby 
galaxy cluster surveys have large datasets of redshift informa- 
tion (see e.g. ENACS Biviano et al. 2004; or 2dFGRS cluster 
sample Burgett et al. 2004) and the redshift information is even 
more limited for galaxy clusters at high redshift. Nevertheless, 
recent surveys of galaxy clusters at high redshift hint that 
the fraction of clusters with substructure is larger at higher 
redshift, thereby suggesting that the hierarchical growth of 
clusters was more intense in the past ( van Dokkum et al. 200C ; 
Huo et al. 200^ 



Haines etal. 2001; 

Halhdav 2005.: .Milvang- Jensen et al. 2008h . 



Demarco et al. 2005 : 



In this paper we have analysed the presence of substruc- 
ture in one of the largest nearby galaxy clusters sample avail- 
able in the literature (Aguerri et al. 2007). Substructure has 
been measured by using the 3D Dressler-Schectman (DS) 
test which, according to numerical simulations, is one of the 
most efficient tests for detecting substructu re in clusters (see 
lKnebe& Miiller 20001: iPinknev et al. 1996 ^. Special attention 
was given to the influence on the results of possible observa- 
tional biases, such as different aperture observations and the 
comparison of different galaxy populations among clusters. We 
were also interested in the study of galaxy substructure in a 
broad range of galaxy densities. Thus, one of the novel re- 
sults presented here is the study of substructure in the outskirts 
(r < 2r2oo) of galaxy clusters, a region that has not been well 
analysed in previous large cluster samples. 

A brief summary of the data is present is Section 2. The 
global substructure is analysed in Section 3. The properties of 
the galaxies located outside and inside substructures are pre- 
sented in Section 4. We present our discussion and conclusions 
in Section 5 and 6, respectively. 

2. Cluster data 

The data used in this study comprises the cluster sample anal- 
ysed in Aguerri et al. (2007; hereafter ASM07). This sample 
consists of 88 nearby (z < 0.1) and isolated galaxy clusters 
with known redshift catalogued by Abell et al. (1989), Zwicky & 
Humason ( 1 96 1 ), Bohringer et al. (2000), or Voges et al. ( 1 999), 
and mapped by the SDSS-DR4 (York et al. 2000; Adelman- 
McCarthy et al. 2006). We only downloaded galaxies located 
within a radius of 4.5 Mpc around the centres of each galaxy 
cluster. Only clusters with a minimum of 30 galaxies with spec- 
troscopic data in the search radius were considered. In order to 
avoid possible incompleteness effects of the SDSS-DR4 spec- 
troscopic data, we completed the spectroscopic SDSS-DR4 ob- 
servations with the data available in the NASA Extragalactic 
Database (NED). Thus, we obtained a constant completeness 
(about 85%) for all magnitudes down to m^ = 17.77 mag (see 
ASM07 for more details). 

Cluster membership in our sampl e was determined us ing a 
combination of the ZHG algorithm ^Zabludoff et al. 19901) and 
the KMM algorithm dAshman et al. 19941) . A first rough cluster 
membership determination was obtained using the ZHG algo- 
rithm. This is a typical gapping procedure that determines cluster 
membership by the exclusion of those galaxies located at more 
than a certain velocity distance (Av) from the nearest galaxy in 
velocity space. We used Av=500 km s"' for our cluster sample 
(see ASM07 for more details). In the second step, cluster mem- 
bership was refined using the KMM algorithm, which estimates 
the statistical significance of bi-modality in a dataset. We ran the 
KMM algorithm on each group of galaxies determined by the 
ZHG algorithm to be galaxy cluster. The KMM algorithm gave 
us the compatibility of the velocity distribution of such group 
of galaxies with a single or multiple Gaussian distributions. We 
considered three different cases: a) Single clusters: the velocity 
distribution of the galaxies is compatible with a single Gaussian; 
b) Cluster with substructure: the velocity distribution show mul- 
tiple groups. We defined the cluster as the most populated group 
plus those groups whose mean velocities lie within 3cr from he 
mean velocity of the largest one; c) clusters with contaminants: 
the velocity distribution is compatible with the presence of sev- 
eral groups, but the mean velocities of the smaller groups deviate 
more than 3cr from the most populated one, which we identify as 
the cluster itself. It is important to note that this cluster member- 
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Table 1. Percentage of clusters with substructure 



Radius 



M,- < -20 Mr < -19 No mag. limit 



r<i'20o 
r200 <r<2r2oo 



11 

55 



33 

57 



34 
44 



Biviano et al. (2002) slightly redefined this test. The param- 
eter 6 was calculated by: 



6i 



1 
a-r 



A 



Nlocal5l 



Si 



(1- ^(Nlocal - 1)/Xt/ 



(2) 



Nhcd 



1)2 



ship selection will aff'ect the detected substructure of these clus- 
ters. We should keep this procedure in mind when comparing 
substructure results from diff'erent studies. The final total number 
of galaxies located in clusters was 10865 objects (see ASM07). 

The 88 galaxy clusters are located in a redshift range be- 
tween 0.02 and 0.1, with an average redshift of 0.07. This red- 
shift range and the fact that the observational data are com- 
plete until nir — \1 .11 mean that the completeness absolute r- 
band magnitude {Mr) is a function of redshift (see Fig. 5 from 
ASM07). Thus, our sample is complete for galaxies brighter than 
Mr - -20.0. It is also important to note that the galaxy data were 
downloaded from SDSS-DR4 archive according to a metric cri- 
terion, and therefore we could be mapping different physical re- 
gions for each cluster. Nevertheless, all the clusters in our sample 
map the region r < 2r2()o (see ASM07), so only galaxies within 
this radius have been considered in the following substructure 
analysis. This provides us with a total number of 6880 galaxies 
located in comparable physical regions for each of our clusters. 

It is important to considered that our cluster membership is 
not free of projection effects. Thus, some galaxies could be lo- 
cated at larger physical radii than the projected ones. Rines et al. 
(2005) showed that ^ 40% of star-forming and 15% of non-star- 
forming galaxies located at projected radii between 1-2^200 were 
located at larger physical distances. 



3. Global substructure in clusters 

The cluster substructure was measured using the kinematic 
and/or spatial information of the galaxies within the clusters. In 
order to account for this, several statistical tests have been devel- 
oped during last decades. In particular, the method developed by 
Dressier & Shectman (1988; hereafter DS) is one of the most ef- 
ficient (Pinkney et al. 1996). The algorithm starts by calculating 
the mean velocity {yi„cai) and standard deviation {criocai) for each 
galaxy of the cluster and its Niocai nearest neighbours. These lo- 
cal values are compared with the mean velocity v^, and standard 
deviation ctc of the global cluster The deviation of the local from 
the global kinematics for each galaxy is then defined by: 



6i 



{ RKVlocal ■ 



cr- 



■ Vcf + (o-h 



'ocal 



(Tcf 



(1) 



Finally, a cumulative quantity A = 25, is computed that 
serves as the statistic for quantifying the substructure. This test 
is normalized with Monte Carlo simulations in which the veloc- 
ities are shuffled among the positions. In this way, an existing 
local correlation between velocities and positions is destroyed. 
The probability of the null hypothesis that there are no such cor- 
relations is given in terms of the fraction of simulated clusters 
for which the cumulative deviation is smaller than the observed 
value. We have normalized the statistic of the test with 1000 
Monte Carlo simulations per cluster. 



with 6v - \viocai - Vc\ and 6^ - maxicTc - o-i„i.ai,Q), where 
the Student-t and x^ distributions were used to calculate the un- 
certainty in the velocity and velocity-dispersion differences, re- 
spectively. This definition of the 6 was designed to obtain groups 
of galaxies that are colder than the cluster and/or have an aver- 
age velocity that differs from the global cluster mean. As be- 
fore, A = 2 6i. This was the version of the test used in our 
study. Nevertheless, we ran the test on our data weighting or not 
weighting 5,. and S^ and obtained similar results. 

The results produced by the test depend on the value of Niocai 
and the number of galaxies per cluster Originally, Dressier & 
Shectman (1988) proposed the computation of A using Niocai = 
10 independently of the number of galaxy cluster members. 
Pinkney et al. (1996) demonstrated that the detection of sub- 
structure was more efficient when Ni„cai depends on the num- 
ber of cluster members. They proposed Ni„cai = ^JNgal, Ngai 
being the number of cluster members. Numerical simulations 
also show that the DS test effectively detects substructures in 
clusters with more than 30 galaxies. Here, we have studied the 
substructure of those galaxy clusters with Ngd > 30 by us- 
ing Niocai - ^JNgal■ Biviano et al. (2002) computed 6/ for each 
galaxy as the average of the (J-values of its N/ocn/-! neighbours. 
We did not find any important differences when carrying out this 
smoothing, so we finally decided to work with the unsmoothed 
5-values of each galaxy. 

3.1. Testing the substructure code 

In order to test the code used for detecting substructure in our 
galaxy clusters we ran several Monte Carlo simulations of clus- 
ters. These simulated clusters did not hold substructure and were 
used in order to test the false positives cases detected by our 
code, i.e. clusters with substructure reported by the test but with- 
out substructure. We considered all our clusters with more than 
30 galaxies in r < 2r2()o, and substituted their velocity distribu- 
tion by a Gaussian distribution with the same Vc and ctc as the 
real clusters. This kind of simulation breaks the possible corre- 
lation between galaxy position and velocity, i.e. possible sub- 
structure. We ran the DS test on these simulated clusters as for 
the real ones, giving no signal of substructure in all cases. This 
means that our test does not detect cases of false positives in the 
study of the global substructure in clusters. 

We have also ran another set of MC simulations of clus- 
ters in order to analyse the efficiency of the DS test for de- 
tecting substructure. In contrast to the previous set of sim- 
ulations, these simulated clusters held substructure. Thus, 
the new simulated clusters had two populations of galax- 
ies: one with no correlation between position and velocities 
showing a Gaussian velocity distribution (we call this a re- 
laxed population) and another of galaxies showing substruc- 
ture (we call this an unrelaxed population). The surface dis- 
tribution of the galaxies of the relaxed and unrelaxed com- 
ponents was simulated following a Hernquist distribution (see 
Hernquist 1990; Mahdavi et al. 1999; Rines & Diaferio 2001^. 
The relaxed galaxies were located within a radius of r/r2oo ^ 2, 
while the unrelaxed component represents a compact galaxy 
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group with radius r/r2oo ^ 0.2, being the centre of the unre- 
laxed component randomly chosen. The velocity distribution of 
the relaxed and unrelaxed components is Gaussian with differ- 
ent mean and sigma. We have selected 7 different values per 
simulated cluster for the relative velocity of cluster and group: 
(vc-Vg)/o-c= 0,0.5, 1.0, 1.5,2.0, 2.5, and 3.0, Vc and 0-^ being the 
mean velocity and velocity dispersion of the relaxed component, 
and Vg the mean velocity of the unrelaxed galaxy population. We 
also varied the ratios between the velocity dispersion of the re- 
laxed and unrelaxed components, choosing o-g/cTc = 0.3, 0.5 and 
0.8 for our simulations. We simulated two sets of clusters with 
different numbers of galaxies (50 and 100). The unrelaxed com- 
ponent in each of these clusters represents 10% and 30% of the 
total galaxy population respectively. Thus, for a fixed total num- 
ber of galaxies we have simulated 42 different clusters varying 
Vg - Vc, cTglcTc and the percentage of substructure. Each of these 
simulated clusters was simulated 100 times changing the cen- 
tre of the unrelaxed component. In total, 8 400 different clusters 
were simulated. 

We computed the fraction of simulated clusters with sub- 
structure provided by the DS test. Figure [1] shows this fraction 
as a function of (vv - Vg)/o-c- Note that the DS test is not very 
efficient for detecting substructure in clusters when the relaxed 
and unrelaxed components have similar mean velocities. In con- 
trast, the efficiency is very high (more than 80%) for clusters 
with (v'c - Vg)/crc > 1-1.5. Note also that, for a given value of 
(vc - Vg)/o-c, the substructure is detected more eficiently when 
it is located in cooler groups, and represents a larger fraction of 
the galaxies in the cluster. These constraints were also pointed 
out by Pinkney et al. (1996) in their numerical simulations. 

3.2. Results on the global substructure 

We applied the DS test to our galaxy cluster sample. We anal- 
ysed the substructure of the clusters, taking into account the 
galaxies inside two apertures: r < r2oo (inner cluster regions) 
and r2oo < r < 2r2()o (outer cluster regions). In the inner cluster 
regions, 34% of our galaxy clusters showed substructure at the 
95% significance level . This value is compatible with other re- 
sults obtained previously using the same statistics. Thus, Solanes 
et al. (1999), Biviano et al. (1997), Escalera et al. (1994), Bird 
(1994), and Biviano et al. (2002) found substructure in 31% 
(21/67), 40%, 38% (6/16), 44% (1 1/25), and 40% (9/23) of their 
clusters respectively. In contrast. Dressier & Shectman (1988) 
found a larger percentage of clusters with substructure: 53% 
(8/15). The amount of substructure in the outskirts of our clus- 
ters is larger. In this case, 44% of the clusters show substructure 
at the 95% significancd level. 



3.3. Substructure and galaxy population 

The value of A was computed taking into account the local en- 
vironment of each galaxy in the clusters. Nevertheless, different 
galaxy populations could show different substructure properties. 
Therefore, one individual cluster could show substructure -or 
not- depending on the galaxy population used to computed A. 
Thus, the study of the substructure in clusters could be biased 
if the clusters in the sample does not have a uniform absolute 
magnitude completeness. 

Our galaxy cluster sample does not show a uniform absolute 
magnitude completeness (see Fig. 5 in ASM07). All the clusters 
have data for the galaxy population brighter than M,- = -20.0, 
but fainter galaxies are lost as the redshift increases. Therefore, 



different clusters trace different galaxy populations. In order to 
analyse the influence of this on the study of the substructure we 
ran the DS test on two galaxy population samples. The first con- 
sisted of those galaxies brighter than M,- - -20.0 located in 
clusters at z < 0.1. The second galaxy population is formed by 
those galaxies brighter than Mr = - 19.0 located in the clusters at 
z < 0.07. In accordance with fig. 5 of ASM07, these restrictions 
allow us to have two families of clusters with uniform galaxy 
populations. 

Table 1 shows the percentage of clusters with substructure 
in their inner and outer regions, taking into account the previous 
restrictions on redshift and galaxy absolute magnitudes. We can 
see that, independently of the galaxy population, the percentage 
of clusters with substructure is higher when the galaxies in their 
outskirts are considered. In the inner cluster regions (r < r2oo) 
few clusters (11%) show substructure if we consider galaxies 
brighter than Mr = -20.0. That percentage increases to 33% 
when a fainter galaxy population is considered. This indicates 
that the detected fraction of clusters with substructure in the in- 
ner regions strongly depends on the limiting magnitude of the 
surveys. Table 1 also shows that the percentage of clusters with 
substructure does not depend strongly on the galaxy population 
in the outer cluster regions. 

3.4. Substructure and global cluster properties 

We have analysed some global properties of the clusters with and 
without substructure, such as the fraction of blue galaxies {ft), 
cluster velocity dispersion (ctc) and the luminosity difference be- 
tween the 2 brightest cluster galaxies (Ami2). For this compar- 
ison we have considered galaxies within an aperture r < r2oo. 
In the previous subsection, we have seen that mapping differ- 
ent galaxy populations in clusters could affect the percentage of 
clusters with substructure. We should therefore take this into ac- 
count and select a cluster sample with a galaxy population free 
of this bias. We have also selected clusters with more than 30 
galaxies. This makes a final number of 28 and 32 clusters when 
we consider galaxies with M, < -20 and M, < -19, respec- 
tively. 

Table 2 shows the mean values of cr^., fi, and Ami2 for clus- 
ters with and without substructure. These values were obtained 
taking into account galaxies in clusters at z < 0.1 and brighter 
than Mr = -20, and those located in clusters at z < 0.07 and 
brighter than M^ = -19. No differences in the clusters prop- 
erties were observed in the two samples of galaxies (see Tab. 
2). Figure |2] shows the cumulative distribution functions of o-^, 
fh and Aoti2 for cluster with and without substructure. This fig- 
ure also shows the cumulative distribution functions of the num- 
ber of galaxy members in clusters with and without substruc- 
ture. Notice that the number of clusters with substructure is very 
small when galaxies brighter than M,. - -20.0 are considered. 
However, this number is larger when we study galaxies brighter 
than Mr - -19.0. In this case, the Kolmogorov-Smirnoff (KS) 
test does not report statistical differences in the cluster properties 
for clusters with and without substructure. Thus, substructure in 
the inner cluster region does not affect ctc fi, or Am 12. 

The above result is not in agreement with Ramella et al. 
(2007). They found a clear difference between the mean value 
of Ami2 for clusters with and without substructure. The work 
by Ramella et al. was based on a sample of 77 nearby clusters 
(0.04< z <0.07) from the WINGS survey (Fasano et al. 2006). 
The substructure was determined using the DEDICA procedure 
(Pisani 1993, 1996). This procedure was aplied to all galaxies in 
the cluster sample brighter than My - -16. The different galaxy 
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Table 2. Properties of clusters with and without substructure 



Mean M, < -20 

Sub No Sub 



M,. < -19 
Sub No Sub 



cr, (km s-') 588.5+199.5 544.8±199.2 567.4±250.7 512.5+150.3 
ft 0.37+0.06 0.42+0.12 0.36+0.10 0.36+0.08 

Amnimag) -0.45+0.5 -0.38+0.4 -0.39+0.47 -0.36+0.34 
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Fig. 1. Fraction of detected clusters with global substructure as a function of (ve-Vg)/(Tj, forMC simulated clusters with (Tc/o-g =0.3 
(left panel), 0.5 (central panel), and 0.7 (right panel). In all panels the symbols represent clusters with a population of re- 
laxed :unrel axed component given by: 90:10 (full points), 70:30 (diamonds), 45:5 (full line) and 35:15 (dashed line). See text for 
more details. 



population studied could be the reason of the disagreement with 
Ramella et al. (2007). 

4. Galaxies in substructure 

Instead of summing all individual values of 6 to get an esti- 
mate of the amount of substructure in the clusters as a whole, 
we can use the values of 6 for each galaxy in order to select 
those galaxies located in substructures. The problem is to deter- 
mine the value 6c, which effectively separates galaxies inside and 
outside substructures. We have taken into account two different 
approaches to this problem. The first approach is to consider dif- 
ferent values of 5, for each cluster In contrast, we can also take 
a unique value of 6g for all clusters. We tested both methods in 
order to find out which one was optimal. 

We selected the values of 6, for each cluster by compar- 
ing the cumulative distribution functions of the values of 6 for 
the galaxies and those from the Monte Carlo simulations of 
each cluster used for normalizing the DS test. These simulations 
broke the correlation between position and velocity of the galax- 
ies. Therefore, they should provide the expected values of 6 of 



galaxies not located in substructures. Three different values were 
selected for each cluster: the maximum value of 6 obtained from 
the MC simulations {Si^max) and the values where the cumulative 
distribution function of 5 values of the simulations was equal to 
0.99 (5,;99) and 0.999 (5;,99.9)- 

The number of galaxies in individual samples is sometimes 
too small for segregation studies. Therefore, the combination of 
data from many clusters will make the statistical analysis more 
reliable. Indeed, the determination of galaxies in substructure by 
adopting a unique value of5g for all clusters is based on the com- 
bination of data in an ensemble cluster. We built an ensemble 
cluster by normalizing the scales and velocities for each galaxy. 
The radial distance of each galaxy to the cluster centre was there- 
fore scaled by the r2oo value for the corresponding cluster, and 
the relative velocity of each cluster galaxy was normalized by 
the velocity dispersion of the cluster. As pointed out by Biviano 
et al. (2002), the ensemble cluster is made on the implicit as- 
sumption that the distributions of galaxy types are similar in the 
individual clusters. Nevertheless, those distributions can be dif- 
ferent for several reasons, such as different observing apertures 
of the clusters and different galaxy incompleteness. Our ensem- 
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ble cluster is free from aperture bias because all our clusters are 
contributing at all radial distances up to Irioo- But we are not 
free from possible biases due to galaxy incompleteness (see Fig. 
5 in ASM07). This bias was avoided by building two ensem- 
ble clusters. The first was built by those galaxies brighter than 
Mr - -20 located in clusters at z < 0.1 (hereafter ECl), and the 
second ensemble cluster was formed with galaxies brighter than 
Mr - - 19.0 in clusters at z < 0.07 (hereafter EC2). In both cases 
only clusters with more than 30 galaxies within 2r2()o contribute 
to the ensemble clusters. Taking into account all previous restric- 
tions, ECl was built with 2593 galaxies from 44 galaxy clusters, 
and EC2 was formed by 2400 galaxies from 34 clusters. 

A global value for 6g was determined by analysing the cu- 
mulative distribution function of the 6 values of the Monte Carlo 
simulations used in the normalization of the DS test. We selected 
two values of 6g for which the cumulative 6 distribution func- 
tion of the MC simulations was equal to 0.95 (dg^gs) and 0.99 
(dg^gg). For our two ensemble clusters these values turned to be 
(5g,95 =1.98 and 2.03 and 6g,gg =2.7 and 3.0 for ECl and EC2, 
respectively. Galaxies with 6 values larger than these global val- 
ues were considered as galaxies in substructures. Note that our 
<5g,95 value is similar to that adopted by Biviano et al. (2002). 

We have used the MC simulated clusters reported in Sec. 3.1 
for determining the best value of c^c for detecting galaxies in sub- 
structures. The substructure in the simulated clusters were anal- 
ysed by the DS test following the same prescription as the real 
ones. In the simulated MC clusters we know whether the galax- 
ies selected as substructure by the DS test are real or spurious 
substructure galaxies. It is expected that the number of galaxies 
in substructure will strongly depend on the adopted value of 6c- 
The best 6c values should maximize the detection of galaxies in 



substructure and minimize the spurious ones. Figure [3]shows the 
number of galaxies detected as substructure versus the number 
of spurious substructure detections for the simulated galaxy clus- 
ters and for the different values of 6c- The results from simulated 
clusters with 50 and 100 galaxies containing 10% or 30% of 
their population in substructures are shown in Fig. [3] Similarly 
to what was reported in Sec. 3.1, the DS test does not detect in- 
dividual galaxies in substructures when the relaxed and the unre- 
laxed component have similar mean velocities. It is also evident 
from Fig. [3]that studies of substructure adopting a unique value 
of 6c for all clusters would be almost dominated by spurious de- 
tections if 6g^g5 was adopted. The spurious detections are lower 
if 6g,gg is used. Studies of substructure using individual values 
of 6c for each cluster will be almost free of spurious substructure 
galaxies by adopting 6ij„a.x- Nevertheless, the fraction of galaxies 
in substructure detected will be very low. The mean fraction of 
galaxies in substructure detected using 6i^gg is 14% higher than 
that adopted with 5,99.9. In contrast, the spurious detections are 
only 6% higher with 6i,gg than with 6i^gg,g- Based on the previous 
considerations, we adopt 5,99 and ^^,99 for the determination of 
the galaxies in substructure. 

4.1. Fraction of Galaxies in substructures 

We computed the fraction of galaxies in substructure as a func- 
tion of their location in the cluster. Figure |4] shows the fraction 
of galaxies located in substructures for ECl and EC2 ensemble 
clusters. In both ensemble clusters, the fraction of galaxies in 
substructure depends on the 6c value adopted. There is a clear 
dependence of the fraction of substructure and the location of 
the galaxies in the cluster. In both ensemble clusters, the fraction 
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of substructure in the inner cluster regions (r < r2oo) is smaller 
than in the regions located at r > r2oo. The EC2 ensemble clus- 
ter shows, independently of the 6c value adopted, a higher frac- 
tion of substructure (especially in the outer cluster region). This 
shows that the fraction of substructure grows as faint galaxies 
are included in the samples. Note that we have not observed any 
segregation between galaxies in substructures and magnitudes 
(see Fig. Is). This means that the fraction of faint galaxies in sub- 
structures is not higher than the fraction of bright ones. 

We have investigated the dependence of the previous results 
with the number of galaxies per cluster In particular, we have 
considered those clusters with more than 50 galaxies and we 
have obtained the fraction of galaxies in substructure as a func- 
tion of distance to the cluster center and galaxy absolute magni- 
tude. No significant diferences were found for this new cluster 
sample. The only difference was a small increase in the fraction 
of galaxies in substructure in the outermost regions of the clus- 
ters (r > r2oo)- The largest variation was for the EC2 sample. 
In this case, the fraction of galaxies in substructures peaks at 
a: 25% in the outermost regions. 

Figure |6] shows the spatial and velocity-radial distributions 
of galaxies inside and outside substructures for the two ensem- 
ble clusters. Figure|6]also shows that the galaxies in substructure 
selected with the 6c values adopted here are mainly located in 
the outer regions of the clusters. We have classified the galax- 
ies into two groups according to their u-r colour. Those galaxies 
showing u-r< 2.22 were called blue galaxies and those with u- 
r> 2.22 were denoted as red. According to Strateva et al. (2001), 
this colour cut separates early- (red) and late-type (blue) galax- 



ies. Note that galaxies in substructures located inside r2oo are 
mainly red galaxies (about 80% and 70% for EC 1 and EC2 re- 
spectively). In contrast, the population of red galaxies in sub- 
structure located at r > r2oo is much lower, but still significant 
(about 60% and 50% for ECl and EC2 respectively). The popu- 
lation of blue galaxies inside raoo is about 20% for ECl and 30% 
for EC2. Nevertheless, the percentage of these blue galaxies lo- 
cated in substructures is very low (less than 5% in all cases). 

The above percentages of blue and red galaxies can be com- 
pared with those obtained from isolated galaxies with similar 
magnitudes. We have measured these percentages using the sam- 
ple of isolated galaxies obtained from Allam et al. (2005). Thus, 
61% and 66% of the isolated galaxies brighter than M,- - -20 
and Mr = -19 turned to be blue ones, respectively. These per- 
centages are larger than those obtained previously for cluster 
galaxies. Notice also than the percentage of blue isolated galax- 
ies is also larger than the percentage of blue cluster galaxies 
located at r > r2oo and in substructures. This could indicate 
that galaxies in substructure are not genuine field galaxies. They 
could be a mixed population of field and cluster galaxies (see the 
discussion about the dynamical state of galaxies in substructure 
in Sec. 5.4). 



5. Discussion 

5. 1 . Spurious Substructure Detections 

In the previous sections we have seen that there is a clear differ- 
ence between the substructure presented in the inner (r < r2oo) 
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and outer (r > r2oo) regions of the clusters, most of the sub- 
structure being located in the outermost regions. It is possible 
that this difference between the inner and outer regions of the 
clusters might not be real, being due instead to our implemen- 
tation of the DS test. We have investigated this effect using the 
MC simulations of clusters shown in section 3.1. In particular, 
we have analysed the fraction of substructure as a function of 
cluster radius detected in the MC simulations of clusters built 
without substructure in them. We need to ensure that the large 



amount of substructure detected in the external regions of the 
clusters (r > r2()o) is not spurious. 

Figure Q shows the fraction of substructure detected in the 
MC-simulated clusters without substructure as function of clus- 
ter radius. The substructure was detected by the DS test as in 
the real clusters using a global value of ^(,,99 or individual values 
for the different clusters 5,-,99 (see section 4). Figure|2]shows that 
the fraction of galaxies in substructure detected at all radial dis- 
tances is between 1 and 2 per cent. There is no excess of galaxies 
detected in substructure in the outer regions of the cluster. This 
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rules out the possibility that the substructure detected in the out- 
ermost regions of the clusters might be dominated by spurious 
detections. 



5.2. The case of Abell 85 

One of our richest clusters is Abell 85, located at z=0.055 with 
273 confirmed members brighter than m,.= 17.88. The substruc- 
ture of Abell 85 has been previously studied in the literature 
(Ramella et al. 2007; Bravo-Alfaro et al. 2009), and there are 
also available X-ray data for this cluster (Durret et al. 2003), 
making this cluster ideal for comparing the substructure obtained 
by us with other studies. 

Bravo-Alfaro et al. (2009) used the DS test to detect the 
substructure in Abell 85 and found five prominent regions of 
substructure. The first was located near the centre of the cluster 
and was identified as C2. Another substructure was found to the 
south called SB, and two more appear in the south-east region 
of Abell 85: one along the X-ray filament detected by Durret et 
al. (2003) was labelled F, and the other along an extension of the 
latter called SE. They also detected a final substructure to the 
west of the cluster (W). Ramella et al. (2007), using a different 
approach, detected three main structures in Abell 85. The two 
substructures (Si and S2) detected by Ramella et al. (2007) are 
close to the C2 and SB substructures reported by Bravo-Alfaro 
et al. (2009). 

Figure [8] shows the substructure in Abell 85 obtained by 
our method. The substructure galaxies were those with 6 > ^g.gg 
and 6 > 61,99 for the M, < -20 (ECl) and M, < -19 (EC'2) 
samples. Note that in the ECl and EC2 samples less substruc- 
ture is measured when 5c = Sg,99- Only in the case of ECl 
are galaxies in substructure obtained in the central region of 
the cluster. These galaxies are part of the C2 and S2 substruc- 
ture detected by Bravo-Alfaro et al. (2009) and Ramella et al. 
(2007), respectively. The substructure measured by this method 
should thus be taken as a lower limit of the real value. The 
situation is better when the galaxies in substructure are deter- 
mined by (5,,99. In this case, for the EC2 sample we obtain two 
groups of galaxies in substructures that correspond to the main 
groups of galaxies in su bstructure (C2 and SE) identified by 
(iBravo-Alfaro et al. 20091) in this cluster. That we have not iden- 
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tified the other groups proposed by these authors might be re- 
lated to the different galaxy population studied. Their observa- 
tions are ~ 2 magnitudes deeper than our data. Notice that we 
have not identified any galaxy from the substructure labeled F in 
Bravo-Alfaro et al. (2009). Indeed, at the position of this group 
of galaxies there are no objects in our data (compare Fig. 4 from 
Bravo-Alfaro et al. 2009 and our Fig. 7). This could indicate that 
the galaxies of this substructure are fainter than Mr = -19 and 
they are not present in our galaxy sample. As we have seen in 
previous sections, the completeness of the observations is im- 
portant for the determination of the number of galaxies in sub- 
structure. Indeed, variations in the abundance of dwarf galaxies 
within substructures have been previously observedas, for exam- 
ple, in the Hercules cluster (Sanchez-Janssen et al. 2005). 

We would Uke also to note that no galaxies in substructure 
were obtained in Abell 85 for the case of EC2 with 6 < ^5,99. 
This might imply the variation of the 6 values from one cluster 
to another Some clusters might have galaxies with larger values 
of 6 that could bias 6g,99 to larger values. This is what happened 
for the EC2 sample; the value of 5^ 99 is so large that clusters like 
Abell 85 do not have galaxies in substructure according to this 
criterion. This means that the galaxies in substructure selected in 
this way would be those with the highest values of 6 and should 
be considered as a lower limit of galaxies in substructures. 

5.3. Luminosity segregation 

The KS test revealed that radial and velocity distributions of 
galaxies brighter and fainter than M^ - -22.0 were statisti- 
cally different (at 99% C.L.). This luminosity segregation was 
also observed in previous studies (see Biviano et al. 2002). In 
most of our clusters, galaxies brighter than Mr - -22.0 are the 
first rank galaxies of the cluster (84% and 59% for ECl and 
EC2 respectively). Those bright galaxies are mostly located out- 
side substructures. Only ^ 3 - 5% (at ECl and EC2) are lo- 
cated in substructures. Moreover, those bright galaxies located 
outside substructures are very clustered and show similar radial 
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velocities to those of the mean cluster (< r/r2oo >~ 0.2 and 
< (v - Vc)l(Jc >~ 0.5 for ECl and EC2). Nevertheless, those 
galaxies with Mr < -22.0 located in substructures are located in 
the outermost regions of the clusters (< r/r20{) >= 1.05 - 1.42 
for ECl and EC2 respectively). 

The previous findings could indicate that those galaxies with 
Mr < -22.0 and not in substructures would be mostly lo- 
cated at the bottom of the cluster potential. In contrast, the 
small fraction of those bright galaxies located in substructures 
would be falling into the clusters associated with galaxy groups 
evolved in a merging process. The origin of these bright galaxies 
{Mr < -22.0) could be due to such accretion or merger processes 
(Governato et al. 2001; de Lucia et al. 2007). 

5.4. Dynamics of galaxies in substructure 

Figure |9] shows the radial and velocity distributions of galaxies 
inside and outside substructures. We have split the galaxy sam- 
ples according to their M-r colours. The red and relaxed galaxies 
are always located closer to the cluster centre than the blue and 
relaxed galaxy sample. Red galaxies outside substructures are al- 
ways cooler than blue and relaxed ones. For galaxies in substruc- 
tures it is also true than the red ones are cooler and are located 
closer to the cluster centre than the blue sample. Nevertheless, 
galaxies outside substructures are basically inside r2oo, in con- 
trast to what happened for galaxies in substructures. 

In earlier sections we saw that, in accordance with our selec- 
tion criteria, our substructure galaxies are mainly located beyond 
r2oo- But can we learn something about the dynamics of these 
galaxies? Are the galaxies in substructure recent aiTivals to the 
cluster potential? Are they an infalling sub-population of galax- 
ies?. These questions can be answered by studying the shape of 
the relative velocity histogram of the galaxies located at r > r2oo- 
Gill et al. (2005) investigated the dynamics of satellite galaxies 
in the outskirts of galaxy clusters from a series of high resolu- 
tion N-body simulations. They found that galaxies in clusters lo- 
cated at r > r20{) were formed by two families: infalling galaxies 
and the so-called back-splash galaxy population. The infalling 
galaxies are entering the cluster potential for the first time. In 
contrast, the back-splash galaxies are located at large cluster dis- 
tances (r > r20{)) but have previously spent time near the cluster 
centre. This back-splash galaxy population could be significant 
in number - up to 50% of the galaxy population located in the 
region 1.4r20{) < r < 2.8r2oo (see Gill et al. 2005). 

The infalling and the back-splash populations can be sepa- 
rated by studying the shape of the relative velocities of the galax- 
ies in the outskirts of the clusters. At large distances from the 
cluster centre (r > 1 .4r2oo) the relative velocity of the infalling 
galaxies is always higher than that of the back-splash galaxies 
(see Gill et al. 2005). Therefore, if the back-splash galaxy popu- 
lation does not exist, then the relative velocity histogram should 
show a Gaussian-shaped peak at relative velocities greater than 
zero. In contrast, the presence of the back-splash population 
should distort the Gaussian shape of the relative velocity his- 
togram, peaking at zero relative velocity (see Gill et al. 2005). 

Figure [TO] shows the relative velocity histograms of all the 
galaxies, and those located inside and outside substructures for 
our ensemble clusters ECl and EC2. We considered only those 
galaxies located at r > 1.4r20{). Figure [TO]indicates that the rela- 
tive velocity histograms of all the galaxies and those located out- 
side substructures peak at zero velocity. Nevertheless, galaxies in 
substructures selected by 5g 99 show a peak in the relative veloc- 
ity histogram different from zero. In order to test the dynamical 
state of the galaxies in substructures, we have compared the rel- 



ative velocities of these galaxies with a mock velocity distribu- 
tion of backsplash plus infalling galaxies from fig. 8 of Gill et al. 
(2005) and Rines et al. (2005). Figure [TT] shows the relative ve- 
locity distribution for a model of backsplash and infalling galax- 
ies. We built several models varying the percentage of infalling 
galaxies. The relative velocity distribution of the models and the 
observed galaxies in substructures were compared using a KS 
test. This provides us the percentage of infalling galaxies located 
in substructure. It should be noticed that the simulations pre- 
sented by Gill et al. (2005) show results for 1 .4r2oo < r < 2.8r2oo- 
Nevertheless, our cluster sample is complete untill 2r2(K). We 
have investigating the variation of the relative velocity distribu- 
tion of galaxies in substructure taking into account those galax- 
ies in the range 2-2.8r2()o for those cluster with galaxies out to 
2.8r2oo- No significant difference was found between the relative 
velocity distributions of galaxies in substructures located in the 
radius ranges 1.4-2.0r2()o and 1.4-2.8r2()o. 

The KS test gives that at the 95% C. L. a pure backsplash 
galaxy population can be ruled out in all cases. Thus, the galax- 
ies in substructure are allways a combination of backplash and 
infalling galaxies. However, galaxies in substructure selected us- 
ing 5,,99 show relative velocity histograms dominated by back- 
splash galaxies. The KS test reports that only 10% for ECl and 
40% for EC2 could be infalling galaxies. In contrast, galaxies in 
substructures selected by <5j_99 are dominated by an infall pop- 
ulation. In this case, the KS test gives that 60% of the galaxies 
in substructures, for ECl and EC2 samples, could be infalling 
galaxies. 

The backsplash scenario in the outskirts of the clusters was 
also studied in the past by Rines et al. (2005) and Pimbblet et al. 
(2006). Similar to the result presented here, they found that the 
outermost galaxy population in clusters is neither a pure back- 
splash population or purely infalling for the first time. They con- 
cluded that local galaxy density is a fundamental parameter for 
the transformation of galaxies in clusters. 

We have analysed the u - r colours of the galaxies in sub- 
structure selected with ^^99 and located at 1.4< r/r2oo < 2.0. 
The blue galaxies (u - r < 2.22) represent 45% and 50% of this 
population of galaxies for ECl and EC2, respectively. As shown 
before, ~ 40% of these galaxies in substructure could be back- 
splash objects. Assuming, that this percentage of objects would 
be red galaxies, then the percentage of blue galaxies of the infall 
populations would be similar to the percentage of blue isolated 
galaxies (see Sec. 4.1). This means that we can not rule out the 
hypothesis that the infalling galaxies would be pure field ones, as 
suggested by recent numerical simulations of cluster formation 
(see Berrier et al. 2009). 

6. Conclusions 

We have analysed the presence of substructure in a sample of 88 
nearby (z < 0.1) and isolated galaxy clusters. The galaxy sample 
consists of 6880 galaxies located in similar physical regions for 
each galaxy clusters (r < 2r2oo)- The substructure was studied 
using the DS statistical test. 

The percentage of clusters with substructure is strongly sen- 
sitive to the galaxy population mapped by the different clusters. 
Thus, 11% and 33% of the clusters show substructure in the in- 
ner regions (r < r2oo) when galaxies brighter than M,- - -20 
or -19 are considered, respectively. The fraction of cluster with 
substructure in their outskirts (r > r2()o) is much larger, be- 
ing 55% and 57% in the two considered cases. No coiTelation 
between substructure and global cluster properties (0-^, fi, and 
Ami2) has been found. 
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Fig. 10. Relative velocity histograms normalized by the velocity dispersion of the clusters for the ECl (top panels) and EC2 (bot- 
tom panels) clusters. The histograms were normalized by their peak, and correspond to all galaxies (left panels), galaxies outside 
substructures (middle panels), and those inside substructures (right panels). All galaxy samples were formed by galaxies located at 
l-4r200 < r < 2.0r2oo- The errors correspond to Poisson statistics. 
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Fig. 9. Spatial (left panels) and velocity-radial (right panels) dis- 
tributions of the galaxies outside (full lines) and inside substruc- 
tures (dotted lines). The galaxies where divided into blue (u- 
r< 2.22; blue lines) and red (u-r> 2.22; red lines). 



We have also studied the galaxies located in substructures by 
the selection of a value of 6c that distinguishes between galax- 
ies inside and outside substructures. The value of 6c was chosen 
by comparing the 5 values measured for the galaxies in the clus- 
ters and those values of 6 from the MC simulations used in the 
normalization of the DS test. This comparison was done individ- 
ually cluster by cluster, or comparing the global distribution of 5 
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Fig. 11. Relative velocity histograms normalized for backsplash 
(dotted line), infalling (dashed line) and total (full line) galaxy 
populations from the Gill et al. (2005) and Rines et al. (2005) 
models of backsplash galaxies. The model showed in this plot 
has 40% of the galaxies as an infalling population. 



values obtained from stacking the data of the individual clusters 
into an ensemble one. In order to avoid possible bias problems 
with different galaxy population traced by the clusters, we have 
also considered two different galaxy populations, one formed by 
those galaxies brighter than M, = -20 located in clusters at 
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z < 0.1 (called ECl) and the other formed by galaxies brighter 
than Mr = -19 in clusters at z < 0.07 (called EC 2). 

The fraction of galaxies in substructure increases with radius 
in both ensemble clusters. Most of the galaxies belonging to sub- 
structures (especially those detected with a global value of 6c) 
are located at r>r2oo- The fraction of galaxies detected in sub- 
structure is higher in EC2 than in ECl. Nevertheless, we have 
seen no trend between the fraction of galaxies in substructure 
and their absolute magnitude. Galaxies brighter than M,. - -22 
are preferentially located outside substructures. Those located in 
substructures (only a few percent) are located in the outer regions 
of the clusters, indicating that they could be located in galaxy 
groups or clusters in the process of merging with the cluster it- 
self. 

We have also investigated the dynamics of the galaxies se- 
lected in substructures. Independent of the method used for se- 
lecting galaxies in substructures, they represent a mixed pop- 
ulation of backsplash and infalling galaxies. The substructure 
galaxies selected by a global value of 6c turned to be domi- 
nated by an infalling population of galaxies. In contrast, those 
galaxies selected in substructures using individial values of 6c 
for each cluster turned to be dominated by back-splash galaxies. 
Assuming that all backsplash galaxies located in substructures 
are red ones, the fraction of blue galaxies of the infalling popula- 
tion is similar to that observed for isolated objects. This indicates 
that we can not ruled out the hypothesis that the infall population 
of galaxies located in substructures would be genuine field ones. 
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